Band Offsets at Semiconductor-Oxide Interfaces from Hybrid Density Functional 

Calculations 
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Band offsets at semiconductor- oxide interfaces are determined through a scheme based on hy- 
brid density functionals, which incorporate a fraction a of Hartree-Fock exchange. For each bulk 
component, the fraction a is tuned to reproduce the experimental band gap, and the conduction 
and valence band edges are then located with respect to a reference level. The lineup of the bulk 
reference levels is determined through an interface calculation, and shown to be almost independent 
of the fraction a. Application of this scheme to the Si-SiCh, SiC-Si02, and Si-HfCh interfaces yields 
excellent agreement with experiment. 
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The discontinuity in the local band structure at 
semiconductor-semiconductor and semiconductor-oxide 
interfaces is a crucial physical property for the operation 
of most electronic and optoelectronic devices [lj. Early 
theoretical research [2J , mainly on semiconductor hetero- 
junctions, provided a deep understanding of theprocesses 
that govern the band alignments at interfaces 0, 0, M, @ • 
This resulted in the development of theoretical models 
providing a reasonable description of band offsets and 
being particularly useful when a wide class of materials 
needs to be screened Q. However, these models mainly 
rely on bulk properties of the interface components [2|, 
and therefore do not account for the detailed atomic and 
electronic properties at the interface, which are known to 
affect band offsets [8j. 

Density functional calculations of band offsets pro- 
vide a qualitative improvement by describing the elec- 
tronic and atomic arrangements at the interface in a self- 
consistent way [H, [l(| U, 12, [HI]. However, the most 
common approximations to the exchange-correlation en- 
ergy, i.e. the generalized-gradient approximation and the 
local density approximation, lead to significant under- 
estimations of band gaps, thereby impairing the relia- 
bility of calculated band offsets. While valence band 
offsets at semiconductor heterojunctions are described 
with reasonable accuracy due to cancellation of errors 
in both interface components, band-offsets errors for 
semiconductor-oxide interfaces can reach several eV jj| 
EH EH, [HI ■ Calculations based on the GW perturbation 
theory yield accurate band offsets at interfaces due to an 
improved description of bulk band gaps [U EH • How- 
ever, these calculations are computationally demanding 
and can only be applied to relatively small systems. For 
instance, the study of realistic semiconductor-oxide in- 
terfaces in which the oxide is noncrystalline [l4, 15, 13] 
is severely hindered. 

In this work, we introduce a scheme for calculating 
band offsets at interfaces through the use of hybrid den- 
sity functionals. These functionals incorporate a fraction 



a of Hartree-Fock exchange [20( and substantially im- 
prove the description of bulk band gaps [H} . We apply 
our scheme to model structures of the Si-SiC>2, Si-HfC>2, 
and SiC-Si02 interfaces, which all feature a realistic de- 
scription of the complex transition region. The band 
structures of the two interface components are lined up 
through their reference levels at the interface [11]. For 
each component, we perform bulk calculations tuning the 
fraction a to reproduce the experimental band gap. The 
lineup of the reference potential in the interface model is 
found to only weakly depend on the fraction a conferring 
consistency on our scheme. For the three interfaces stud- 
ied, the calculated band offsets are in excellent agreement 
with experiment. 

We considered a class of hybrid density function- 
als based on the generalized gradient approximation of 
Perdew, Burke, and Ernzerhof (PBE) [22], which are ob- 
tained by replacing a fraction a of PBE exchange with 
Hartree-Fock exchange [23[. The functional defined by 
a=0.25 is referred to as PBEO and is supported by the- 
oretical considerations *23J. Core- valence interactions 
were described through normconserving pscudopotentials 
generated at the PBE level. The valence wave func- 
tions were expanded in a plane-wave basis set defined 
by an energy cutoff of 70 Ry. The interface calculations 
corresponding to large supercells were performed with 
a Brillouin-zone sampling restricted to the T point. In 
the bulk calculations, the positions of the band extrema 
were determined through converged fc-point samplings. 
The integrable divergence of the Hartree-Fock exchange 
term was explicitly treated [3]. Structural relaxations 
were carried out at the PBE level. We used the imple- 
mentations in Quantum-ESPRESSO 25| and CPMD [26|. 



The choice of the model structures requires particular 
attention. Indeed, previous density-functional studies on 
crystalline-crystalline Si-ZrC>2 and SiC>2-Hf02 interfaces 
have revealed a strong sensitivity of the band offsets on 
the adopted model of the interfacial bonding pattern [§] . 
Therefore, we considered model interfaces in which the 
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TABLE I: Band gaps (in eV) of Si, SiC, Hf0 2 , and Si0 2 , 
calculated using functionals incorporating a varying fraction a 
of Hartree-Fock exchange: a=0 (PBE), a=0.25 (PBEO), and 
an optimal fraction ao reproducing the experimental band 
gap. 





PBE 


PBEO 


Optimal (qq) 


Expt. 


Si 


0.6 


1.8 


1.1 (0.11) 


1.1 


SiC 


2.2 


3.9 


3.3 (0.15) 


3.3 


HfC-2 


4.3 


6.7 


5.9 (0.15) 


5.9 


SiQ 2 


5.4 


7.9 


8.9 (0.34) 


8.9 



oxide is amorphous to ensure that the transition region 
is smooth in terms of bond parameters and coordination 



For the Si-Si02 interface, such a choice led to 
variations of 0.1 to 0.2 eV in the band offsets calculated 
at the PBE level (27|. The Si-Si0 2 interface adopted 
here is described through a 217-atom superlattice model 
in which layers of crystalline Si (9 monolayers) and of 
amorphous SiC>2 (17 A) alternate This interface 

model structure incorporates a set of atomic-scale fea- 
tures inferred from experimental data. The SiC-Si02 in- 
terface is modeled by a 237-atom slab characterized by a 
chemically abrupt transition between crystalline 4_ff-SiC 
(8 monolayers) and amorphous SiC>2 (16 A) 15]. The 
Si-HfC>2 interface is described by a superlattice model 
comprising 282 atoms and including a Si02 interlayer (7 
A) between crystalline silicon (10 monolayers) and amor- 
phous Hf0 2 (12 A) 0. To represent the bulk of the 
interface components, we used the corresponding crys- 
talline structures for Si and 47J-SiC, a disordered model 
for SiC>2 [28| , and the monoclinic structure for HfC>2 [2!| . 
In Table HI the band gaps of these reference bulk com- 
ponents are calculated at the PBE and PBEO levels and 
compared to experimental values. As well known, the 
PBE band gaps severely underestimate their measured 
counterparts. While the inclusion of Hartree-Fock ex- 
change always enhances the calculated band gap, the 
comparison with experiment is not systematically im- 
proved at the PBEO level. 

First, we calculated band offsets at the three interfaces 
within both PBE and PBEO. The bulk band edges were 
aligned through a reference potential calculated across 
the interface [11(. As reference potential, we generally 
used the local potential [3(J, but for SiC>2 and HfC>2 
we resorted to the energy levels of the deep O 2s and 
Hf 5s states which are less sensitive to structural disor- 
der. Focusing on the Si-Si02 interface, we show in Fig. [1] 
the planar-averaged electron density and local potential 
across the interface. With respect to PBE, PBEO only 
yields a small redistribution of the electron density [Fig. 
Gib)], which results in a difference of AV=0.16 eV be- 
tween the lineups of the potentials at the interface [Fig. 
Hid)]. This indicates that the dipole contribution to the 
band offsets is already well described at the PBE level. 
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FIG. 1: Planar-averaged (a) electron density and (c) local po- 
tential across the Si-Si0 2 interface calculated in PBE (solid) 
and PBEO (dashed). The difference between the electron den- 
sities and the local potentials in the two schemes is shown in 
(b) and (d), respectively. The shaded areas correspond to the 
transition regions between Si and Si0 2 . 



TABLE II: Valence (AE V ) and conduction (AE C ) band offsets 
at the Si-Si0 2 , SiC-Si0 2 , and Si-Hf0 2 interfaces calculated 
in PBE, PBEO, and the mixed scheme. Experimental band 
offsets are from Refs. [3lJ, [34]], and [35|, respectively. 







PBE 


PBEO 


Mixed 


Expt. 


Si/Si0 2 


AE V 


2.5 


3.3 


4.4 


4.4 




AE C 


2.3 


2.7 


3.4 


3.4 


4tf-SiC/Si0 2 


AE V 


1.4 


2.0 


3.0 


2.9 




AE C 


1.7 


2.0 


2.6 


2.7 


Si/HfQ 2 


AE V 


2.3 


3.1 


2.9 


2.9 




AE C 


1.5 


1.9 


1.7 


1.7 



Nevertheless, band offsets calculated at the PBEO level 
noticeably improve upon the PBE ones. For instance 
the valence band offset goes from 2.5 eV to 3.3 eV, to 



be compared with the experimental value of 4.4 eV 31 1. 
This improvement is mostly due to a better description of 
bulk band gaps in PBEO (Table [TJ) . However, deviations 
with respect to experiment are still remarkable. Similar 
observations also hold for the other interfaces (Table HI)) . 

To further improve band offsets, it appears impera- 
tive to more accurately describe the band gaps of bulk 
components. This cannot be achieved for both interface 
components through the use of a hybrid functional with 
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FIG. 2: (color online) Valence band maximum and conduc- 
tion band minimum vs a for four different materials. Shifts 
calculated with GW and quasiparticle self-consistent GW 
(QSGW) [l8| are also shown. Vertical lines represent the 
experimental band gaps and are shown in correspondence of 
ao (see text). 




FIG. 3: (Color online) Schematic illustration of band offset 
determination at the interface between components A and 
B through the mixed scheme (see text). The local poten- 
tial V(z) is obtained from an interface calculation based on 
a hybrid functional with a — (a^ + ajf)/2 and determines 
the lineup of the reference levels in the two components. For 
each component, the band extrema are then aligned to the 
reference potential through bulk calculations based on hybrid 
functionals with a material-specific ao chosen to reproduce 
its experimental band gap. AE V and AE C are the resulting 
valence and conduction band offsets. 



a fixed fraction of Hartree-Fock exchange a. However, 
it has been argued that there is no universal fraction a 
for all materials and that its optimal value could even 
be property dependent [3^ ]. In Fig. [21 we show the evo- 
lution of the conduction and valence band edges of the 
four bulk materials considered in this work as a function 
of a. Since the bulk band gaps increase linearly with a, 
the value of a can always be chosen to reproduce the ex- 
perimental band gap. This resulted in optimal values ao 
of 0.11 for Si, 0.15 for SiC and Hf0 2 , and 0.34 for Si0 2 
(Tableland Fig. [3]). We note that for these a values, 
the hybrid functional calculation not only gives the ex- 
perimental band gap but also the position of the band 
extrema with respect to the adopted reference potential. 

The consideration of a fraction a of the Hartree-Fock 
exchange is equivalent to an effective static screening of 
the long-range interaction: a ~ I/cqo, where eoo is the 
electronic part of the dielectric constant. Indeed, the op- 
timal ao's found above respect this relationship in a qual- 
itative way. To further support the proposed adjustment 
of a, we compared the band shifts of Si and SiC"2 calcu- 
lated with the hybrid functionals to those obtained with 
GW and quasiparticle selfconsistent GW (QSGW) [Hj]. 
The GW band shifts are reported in Fig. [2] in correspon- 
dence of the value of a which reproduces the band gap 
found in the GW calculation. The comparison between 
hybrid and GW band shifts is quite good, with differences 
not exceeding those between different GW schemes. 

The observations above can be combined to obtain 
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FIG. 4: (Color online) Lineup of the local potential (filled 
circles), expressed as a difference with respect to the result 
obtained with PBE (a—0), vs fraction of Hartree-Fock for the 
Si-Si02, SiC-Si02, and Si-HfO^ interfaces. Squares are shown 
in correspondence of the ao's of the individual bulk compo- 
nents (0.11 for Si, 0.15 for SiC and Hf0 2 , and 0.34 for Si0 2 ). 
Filled symbols correspond to actual calculations, while open 
symbols are positioned by interpolation or extrapolation. 



band offsets. From the interface model, the lineup of the 
reference potentials in the two bulk components is deter- 
mined. Then, bulk extrema are positioned on each side 
of the interface using the results of hybrid functional cal- 
culations with optimal ao for each bulk component [33J|. 
This mixed scheme is graphically illustrated in Fig. [31 

The scheme critically relies on the fact that the lineup 
extracted from the interface calculation is not signifi- 
cantly dependent on the adopted fraction of Hartree-Fock 
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exchange a. To further support this point, we performed 
several hybrid functional calculations with varying a for 
each of the model interfaces under consideration. As 
shown in Fig. [U the lineup between the reference lev- 
els in the two bulk components only marginally depends 
on a, generalizing the observation made above for the 
Si-Si02 interface. The dependence on a is even weaker 
for the other interfaces. Inspection of Fig. QJb) shows 
that the charge transfer occurring at the interface is sim- 
ilar to that observed in the oxide, suggesting that the 
change of the interfacial dipole should be attributed to 
a modification of the local chemistry rather than to the 
increase of the number semiconductor-induced gap states 
Q resulting from the modified band offsets. Assuming 
that the interface calculation is performed with a value 
of a corresponding to the average of the ao's pertaining 
to the two interface components, we estimate that the 
induced indetermination is smaller than 0.15 eV in the 
worst case. 

The band offsets obtained through the application of 
the mixed scheme are included in Table [ill The calcu- 
lated values agree extremely well with experimental ones. 
The error between theory and experiment is similar to 
the intrinsic indetermination of our approach and to the 
scatter between different experiments. In comparison to 
hybrid calculations with a fixed fraction a, the mixed 
scheme provides a striking improvement in the theoreti- 
cal estimation of band offsets. 

In conclusion, we demonstrated the accurate determi- 
nation of band offsets through the use of hybrid den- 
sity functionals and experimental band gaps. The pre- 
sented scheme constitutes a predictive tool which is com- 
putationally less demanding than GW calculations, yet 
achieves band offsets of comparable accuracy [18( . Its ap- 
plication to complex interface components such as non- 
crystalline oxides is within reach without loss of accuracy. 
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